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Mixed three-point correlation functions of the nonlinear integrated 
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ABSTRACT 

In this paper I investigate the family of mixed three-point correlation functions <T 9 y 3 " 9 ), 
q = 0, 1,2, between the integrated Sachs-Wolfe (iSW) temperature perturbation t and the 
galaxy overdensity y as a tool for detecting the gravitational interaction of cosmic microwave 
background (CMB) photons with the potentials of nonlinearly evolving cosmological struc- 
tures. Both the iSW-effect as well as the galaxy overdensity are derived in hyper-extended 
perturbation theory to second order and I emphasise the different parameter sensitivities of 
the linear and non-linear iSW-effect. I examine the configuration dependence of the relevant 
bispectra, quantify their sensitivities and discuss their degeneracies with respect to the cos- 
mological parameters Q m , cr 8 , h and the dark energy equation of state parameter w. I give 
detection sigificances for combining PLANCK CMB data and the DUNE galaxy sample, by 
using a quadratic approximation for the likelihood with ACDM as the fiducial cosmology: 
The combination of PLANCK with DUNE should be able to reach a cumulative signal to 
noise ratio of 0.6cr for the bispectrum (ry 2 ) up to £ — 2000, where the most important noise 
source are the primary CMB fluctuations. 

Key words: cosmology: cosmic microwave background, large-scale structure, methods: ana- 
lytical 



1 INTRODUCTION 

The integrated Sachs -Wolfe (iSW) effect dSachs & Woifelll967l ; 
Hu & Sugiyama 1994) is a secondary anisotropy in the cosmic mi- 
crowave background (CMB) sky due to gravitational interaction 
of CMB photons with time- varying potentials in the cosmic large- 
scale structure (LSS) . The linear iSW-effect is an interesting obser- 
vational channel because it directly measures the dark energy (DE) 
content of the universe due to its influence on the growth function. 
The linear iSW-effect has been det ected with high statistical signif- 
icance with different tracer obje cts (Nolta et al. 2003; F osalba et al.l 
bOOSl fciannanto nio et alj20 06: Rassat et al. 2007), and constraints 
on cosmological parameters such as Q,„ and the dark energy equa- 
tion of state parameter w can be derived ( Crittenden & Turo kl 19961 : 
iGiannantonio et al.ll2008l : lDouspis et~alll2008h . 

On smaller angular scales one expects a contribution to the 
iSW-effect caused by gravitational interaction of the CMB photons 
with nonlinearly evolving cos mic structures, which i s refered to as 
the Rees-Sciama (RS) effect dRees&Sciarnll 19681) . The angular 
power spectrum of the RS-effect has b e en derived in perturbation 
theory ^Martinez-Gonzalez et alJll99Ct iTuluie et al.lll996l : fSeliakl 
ll996l : ICooray|2002l : ISchafer & BartelmannkOOfj). and model led on 
n-body simulations (Puchades e t al.ll2006l : lMaturi et alj2007l) . 



In this paper, I intend to derive mixed bispectra between the 
iSW-temperature perturbation and the density of tracer objects, as 
a tool of measuring the RS-effect caused by nonlinearly evolving 
cosmic structures, which give rise to non-Gaussian features in the 
CMB. I use hyper-extended perturbation theory for describing non- 
linear growth and summarise the key equations related to cosmic 
structure formation in Sect. [2j In addition, I discuss the different 
parameter sensitivity of the RS-effect compared to the iSW-effect 
in this section. The mixed bispectra (y 3 ), (y 2 r) and (yr 1 ) of the 
galaxy overdensity y and the iSW-temperature perturbation t are 
derived in Sect, fj] and I provide a discussion of their configu- 
ration dependence. I focus on these three bispectra because (t 3 ) 
has already been shown to be very small ( Mollerac h et al.l 1 19951 : 
ISpergel & Goldberdll999l : IVerde & Spergelll2002h . and because it 
would not be possible to distinguish primordial non-Gaussian fea- 
tures in the CMB from those induced by a secondary anisotropy 
without adding information about the LSS in form of the tracer 
density. In Sect. [4] I derive signal to noise ratios for the measure- 
ment of these bispectra for combining PLANCK^] CMB data with 
the main galaxy sample of the Dark UNiverse Explore^] (DUNE), 
which covers half of the sky and contains about 3 x 10 9 objects. 
Then I continue with a discussion of the parameter sensitivity of 
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the mixed bispectra (the relevant parameters being Q m , <r 8 , w and 
h) and summarise my main results in Sect. [5] 

At this point I should emphasise that I consider only idealised 
measurements, i.e. the PLANCK CMB -observation is tainted with 
a Gaussian noise component and a Gaussian beam, and the galaxy 
surveys are assumed to have a Poissonian noise component in the 
galaxy number density. A quantification of systematic effects, such 
as unresolved point sources or the kinetic Sunyaev-Zel'dovich ef- 
fect, which both are equally associated with overdensities in the 
LSS, is beyond the scope of this paper. 

The cosmological model used is the spatially flat ACDM 
cosmology with adiabatic initial conditions. Parameter values are 
Q.„, = 0.25, <r 8 = 0.8 and n„ = 1, with a constant equation of 
state w = — 1 for the dark energy fluid. The Hubble-constant has 
the value Ho = 100 h km/s/Mpc. For simplicity, the DUNE galaxy 
sample is assumed to have a non-evolving bias of unity, b = 1, 



2 KEY FORMULAE 

2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the matter density 
parameter Sl m , the Hubble function H(a) = dlna/df is given by 



H 2 (a) n„, I r l 

-^ + (l-n,„)expl3 J dlna(l + w(a))\, ill 



Hi 



with the dark energy equation of state w(a). The value w = -1 
corresponds to the cosmological constant A. The relation between 
comoving distance x an d scale factor a is given by 



X = c 



J da — , 
a 2 H(a) 



(2) 



with the speed of light c. 



2.2 Linear structure formation and power spectra 

For the linear CDM density power spectrum P(k), defined from the 
fluctuation amplitude of the density field, (d(k)d(k')) = (2n) 3 S D (k+ 
k')P{k), I use the the ansatz 



P(k) cc k"'T 2 (k), 



(3) 



with the transfer function T(k). The transf er function is appro xi- 
mated with the polynomial fit proposed by ( Bardeen et al. 1986), 



T(q) 



ln(l + 2.34?) 
23Aq 



(l + 3.89^ + (\6.lq) 2 + (5A6q) 3 + (6.71#) 



The wave vector k = qT is rescaled with the shape parameter T 
JSugivamalll995h . 



r = Q m h exp -Cl b 1 + 



V2A 



(4) 



The spectrum P(k) is normalised to the variance cr 8 on the scale 
R = 8 Mpc/h, 



t2r = dk^Pi^W'ikR), 



(5) 



with a Fourier transformed spherical top hat filter function, W(x) = 
3 ji(x)/ x. j((x) is the spherical Bess el function of the first kind of 
order t ( Abramowitz & Stegun 1972). The homogeneous growth of 
the density field in the linear regime, 6(x, a) = D+{a)6(x, a = 1), is 



described with the grow th function D+(a), which results from solv- 
ing the growth equation l Wang & Stein hardt 1998; lTurner & White! 



ll997l ; lLinder & Jenkinsll2003l) . 



+ - (3 + ^-^) ^-D + (a) = ^n m (a)D + (a), (6) 
act 1 a\ am a J aa 2a l 

and which assumes the simple solution D + (a) = a in the SCDM 
cosmology, where 0,„ = 1 and H{a) oc cr 3/2 . 



2.3 Galaxy biasing 

The fractional perturbation An/(n) of the mean number density {n} 
of galaxies is related to the overdensity 6 of dark matter. I use the 
phenomenological relation 



An 
(n) 



Ap 
'<P> 



■■ bd, 



(7) 



with a constant bias parameter b. There exist more elaborate non- 
linear biasing models, even including time evolution (which itself 
can be parameterised with e.g. b(a) = b + (1 - a)b a ), but for sim- 
plicity, I will work with a constant linear biasing model. 



2.4 Linear and nonlinear iSW-effects 

The iSW-effect is caused by gravitational interaction of a CMB 
photon with a time-evolving potential O. The frac tional perturba- 
tion t of the CMB temperature Tcmb is given by dSachs & Wolfd 
[1963) 



AT 



c 1 J dq c 5 Jo da 



(8) 



where q denotes the conformal time, and the integration is extended 
to the horizon distance^. In the last step, I have replaced the inte- 
gration variable by the comoving distance^-, which is related to the 
conformal time by Ax = -cAq = -cAt/a, and the time derivative 
of the growth function has been rewritten in terms of the scale fac- 
tor a, using the definition of the Hubble function Aa/At = aH(a), 
with the cosmic time t. The gravitational potential O follows from 
the Poisson equation in the comoving frame, where the Newton's 
constant G is replaced with the critical density p cl j t = 3H 2 /(8nG), 



AO : 



3H 2 Q„, 
2a 



(9) 



Substitution yields a line of sight expression for the linear iSW- 
effect t {1) (integrating along a straight line and using the flat-sky 
approximation), sourced by the linear density field <5 (1) , 



3£X 



C XH , A (D + \ A" 1 m 

1 ^^(t)^' 



(10) 



with the inverse (dimensionless) Laplace operator A 1 ld 2 H solving 
for the (dimensionless) potential if, 



"h 



(11) 



with the Hubble distance d H = c/H , Extending this expression 
to include a nonlinear correction 6 {2) to the density field, 6 = 
D + (a)d m + D 2 + (a)S< 2 \ gives a contribution t (2) from nonlinear struc- 
ture formation, 



r (2) 



A X a l H(a) — — -rr<5 • 

c J D Aa \ a I d- H 



(12) 
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Figure 1. Time evolution Q y (a) and Q T (a) of the source fields as a function 
of scale factor a: density field in first order tfW (solid line), density field in 
second order <S' 2 ' (dashed line), potential derivative in first order </>''* (dash- 
dotted line) and potential derivative in second order </>' 2 ' (dotted line). 



This second order effect dRees & Sciama|[l968l) has in fact differ- 
ent parameter dependences, as illustrated by a simple example: 
The iSW-effect vanishes in the SCDM-model (fl„, = 1), because 
D + (a) = a, and is nonzero in models with dark energy, because the 
growth function D+(a) scales slower at low redshifts. On the con- 
trary, the RS-effect is present in the SCDM-cosmology, as it mea- 
sures d(D+/a)/da instead of d(D + /a)/da, and is at low redshifts 
always smaller in dark energy cosmologies compared to SCDM . 

The line of sight expression for the projected galaxy overden- 
sity y in first and second order is given by 

dz 



X 

Jo 



d X p(z)^D + b6 w , 
d X p(z)—D 2 + b6 <2 \ 



(13) 



(14) 



with the redshift distributio n p(z)dz and the li near bias parameter b. 
p(z)dz is approximated by dSmail et al. 1995), 



p(z)dz = po I — | exp 



-l-l Idz 

zo. 



with l = !£r(i 

Po P \P 



(15) 



Fig.Q]shows the time-evolution Q y (a) of the the source fields S {1) . 



Q y {a) 



D + (a), 
D 2 + (a), 



for<5 (1) 
for,?* 2 ' 



(16) 



and Q T (a) of y> (,) which result in the galaxy overdensity y and for 
the iSW-effect r by projection, respectively: 



(17) 



I define Q q (d) to be the time evolution of the bispectrum 
{<p q 8^~ q } of the source fields, which acquires terms from first and 
second order growth (anticipating results from Sect. l3.3t . 

3Dt for<<5 3 ) 

for (6 2 (p) 



Q q {d) 



^(i D i)(£^)^l(^)\ for^ 2 > 
3(i^) 2 (£#) for<^). 



(18) 



Figure 2. Time evolution Q q {a) of the bispectra (r''y 3_<? >, for (y 3 ) (solid 
line), (ry 2 ) (dashed line) (r 2 y) (dash-dotted line) and (r 3 ) (dotted line), in 
second order perturbation theory, as a function of scale factor a. 



It is apparent how the different terms in the first and second order 
time evolution of the density and the potential field affect the time 
evolution of the bispectrum {tp q 6 3 ~ c '). Up to this moment I have 
included the iSW-bispectrum S3 = (t 3 ) for completeness. I would 
like to point out that it is not possible to distinguish non-Gaussian 
features imprinted into the CMB by the iSW-effect from primordial 
non-Gaussianities, due to the achromaticity of the iSW-effect. For 
this reason, I will only consider the bispectra (y 3 ), (ry 2 ) and (T 2 y), 
which can be measured in cross-correlation with a tracer field 7. 
Fig. |2 shows the time-evolution Q q (a) of the 3-point correlation 
functions {tp q 5 i ~ q } as a function of scale factor a. 



2.5 Weighting functions 

For composing the weighting functions W q (x) needed in the pro- 
jection l |35l l, one can read off the weightings 



(19) 



(20) 



W T (x) = ^-a 2 H(a), 
c 

W r (j() = bp(z)^ = -p(z)H(a) 
d X c 

from eqns. ([8} and d!4t . which both have units of inverse (Mpc/h). 
With W T (x) and W y (x), the weighting functions W q (x) used for car- 
rying out the projection of the soure field bispectra (tp q S 3 ~ q ) to the 
angular bispectra {r q y i ~ q ) can be written down, 

for <y 3 > 



Wfy) = W q (x)Wl->(x) ■■ 



for (y 2 T) 
for <yr 2 > 
for <r 3 >. 



(21) 



W 2 (x)W T (x), 
W y (x)W 2 (x), 

The weighting functions W^Ct) are depicted in Fig. [3] Common to 
all W q (x) is a broad peak between 1 Gpc//z and 4 Gpc/fr, corre- 
sponding to the maximum of the galaxy redshift distribution. 



3 MIXED 3-POINT FUNCTIONS 
3.1 Bispectra 

Bispectra are a tool for quantifying non-Gaussianities which arise 
in the course of structure formation when nonlinearities in the 
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Figure 3. Line-of-sight weighting functions W q (x) for the projection of the 
bispectra {r q y i ^ q ), for (y 3 ) (solid line), (ry 2 ) (dashed line) and (r 2 y) (dash- 
dotted line), as a function of comoving distance^, for the redshift distribu- 
tion of the DUNE galaxy sample. 



structure formation equations set in. Bispectra are the Fourier ana- 
logue of three-point correlation functions, i.e. for the density field 
8 the relation 



(6(kOS(k 2 )S(k 3 )) = (2n) 3 d D (k l + k 2 + k 3 )B 6 (k u k 2 , Jfc 3 ) 



(22) 



defines the bispectrum S, s via the three-point variance of the 
Fourier-transformed density field. Equivalent formulae apply to an- 
gular bispectrum B y of the projected density field y 

WWaM«) = Qx?s D (fi + e 2 + e 3 )B r ({ 1 ,{ 1 , e 3 ), (23) 

in the flat-sky approximation, with angular wave vectors t. 



3.2 Notation 

I introduce a notation reminiscent of weak lensing tomography: 
The bispectra B( a &c)(A,&>^3) are indexed by the variables a,b,c, 
which assume the value one for the iSW-field r and zero for the 
galaxy overdensity y. This is a way to distinguish between e.g. 
the bispectra B m 0u ( 2 ,h) = (rtfiM^Ws)), B m) (.t u { 2 ,£ 3 ) = 
(r(AM4M4)> and B m) {€ u € 2 ,€3) = <r(^)r(4M^)>, assuming 
Yt P ( P = 0. These bisepctra are not equal (meaning the triangles can 
not be mapped onto each other for general configurations), unlike 
permutations in weak lensing tomography bin numbers. Most of the 
signal in an actual measurement, however, will come from general 
configurations and not from isosceles or equilateral configurations, 
for which these bispectra would be in fact identical. Nevertheless, 
all permutations with fixed q = a + b + c share the same time evo- 
lution Q q (a) and the same weighting function W q (x). 

I relate the source field bispectra B (abc) to the bispectrum of 
the density field S ( ooo) with the formula 



B (ahc) (ki,k 2 ,k 3 ) = I — 
\<F H 



1 



K l K 2 *3 



B(ooo)(k],k 2 , k 3 ), 



(24) 



and average over all combinations with fixed q for visualisation 
purposes: The bispectra B q {t\,t 2 ,{ 3 ) depicted in Figs. [4] through [7] 
will be averages over all permutations, i.e. 

(-D ? 



Bq(tl,l 2 ,( 3 ) 







B(abc)(fl>^2> ^3)1 



(25) 



with the binomial coefficient ( 3 J = -^zm counting the permuta- 
tions with a fixed number of q iSW-fields t and 3 — q galaxy over- 
density fields 7. In total, there are three permutations for B\ and 
B 2 , and one permutation for So and S3. The occurence of first and 
second order terms of the respective fields S, <p in the perturbative 
expansion determines the time evolution of the bispectrum (ip^S 3 ^' 1 ) 
and lead to the equations in Sect. 12.41 



3.3 Nonlinear structure formation and perturbation theory 

In order to carry out the evaluation of the bispectra in perturba- 
tion theory I use the reformulation of the iSW-effect outlined above 
such that the souce fields and are dimensionless, as well as 
their respective first and second order time evolution Q T (a), Q y {a) 
and their second order perturbations S (2) and tp {2) . The weighting 
functions W r 0f) and W y ix) used for carrying out the line of sight 
projection will have units of inverse Mpc//j, such that the units 
are consistent in the Limber-projection of bispectra, and that the 
resulting ang ular bispectra {T q y i ^ q } are dimensionless. Using the 
formalism of Seliak ( 1996), I separate the time evolution and the 
mode-coupling of the source fields. 

The first order contribution to the bispectrum B (k[ , k 2 , k 3 ) = 
((5(Ai)(5(A:2)'5(^3)) of the density field from nonlinear structure for- 
mation is can be expanded into a pr oduct of power spectra Pss(k) 
(I will follow iTakada & Jainll2003aL and use th e nonlinear power 
spectrum for Psg(k)) and is given by jFrvlll984a llbl): 



B (k l ,k 2 ,k 3 )= Y> M(ki,kj)P S6 (ki)Pss(.kjX 



(26) 



a.m 1,2.31 

i*3 



with the mode coupling functions of classical perturbation theory, 

=t + (HM* (27) 

where x = kjkj/ (kjkA denotes the cosine of the an- 
le between k t and In hyper-extended perturbation theor y 



Scoccimarro & Frieman|[l999l : IScoccimarro & Couchmanl l200lt) . 
the mode coupling function is replaced by: 

M(ki,kj) = ^aikdaikjHbik.Mkj)^ + ^)x+lc(.k i )c(k J )*?.(2$) 



The coefficients a(k), b{k) and c(k) are given by 
a(k) 



1 + o-- a2 (z) V0.7g(rc) (q/4)" +3 - 5 



b(k) 



c(k) 



1 + (q/4) n+3S 
1 + 0.4(n + 3)q" +3 



1 +q" 



1 + 



4.5 



1.5+(b+3)' 



- Q4f 



1 + (2q) n 



(29) 
(30) 

(31) 



where the time evolution of the fluctuation amplitude is given by 
the linear growth function, crg(z) = D + (z)(Tg. In eqns. l |29l l, d30t 
and J3U . the wave vectors k are expressed in units of the nonlinear 
wave number fc NL , q = k/k NL . The nonlinear wave number at scale 
factor a is given by the scale at which the variance a 2 of the density 
fluctuations becomes unity, 



& 3 kD\ 





(z)P(k) = 1 



Ankl L Dl(z)P(k HL ) = 1. (32) 



The logarithmic slope of the linear power spectrum, 
din S(/t) 



n(k) 



dlnfc 



(33) 
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can be directly derived with form of the transfer function T(k) given 
in eqn. ( 12.2b . and is used for determining the saturation parameter 
Q(n), which is defined as the logarithmic slope n of the linear CDM 
spectrum, 

4-2" 



Q(n) 



(34) 



1 + 2" +1 

I use the functional form derived bv lSmith et alU 2003) for the non- 
linear CDM spectrum Pssik) and its slow time evolution, parame- 
terised with Q.„,(a). 

3.4 Limber projection 

In order to relate the source field bispectra ((p q S 3 ~ q ) to the angular 
bispectra (T q y 3 ^ q ) of the obser vables r and y, I carry out a Limber- 
style projection dLimbedll954l ). 



Jo X 



W q (x)Q q (x) *,(*!, fe,*3), 



(35) 



for which I adopt the flat-sky approximation, which is justified as 
the nonlinear iSW-effect is a small-scale phenomenon. The relation 
between the wave vectors k p and the multipole vectors t p is given 
by k p = t p lx, P = 1,2, 3. The spherical bispectru m B a (l\,l2,l-j) is 
related to the flat- sky bispectrum B q {t\ , f 2 , ^3) via dMiralda-Escudd 
ll99ll : lKaiserlll992h 



(2 




£3 




uU(2f p + 1) 



An 



-BJt u t 2 ,t 3 \06) 



where 

(« 




('2 








dx P h (x)P e2 (x)P (2 (x), 



(37) 



x = cos 6, denotes the Wigner-3 j symbol, which re- 
sults from integrating over three Legendre polynomials Pf(x) 
dAbramowitz & Stegunlll972h . The Wigner-3 j cancels configura- 
tions which do not satisfy the triangle inequality \(i ~ tj\ ^ 4 ^ 
Cj + £/. The factorials arising in the evaluation of the Wigner-3 j 
symbol are computed using the Stirling-approximation for the T- 
function, T(n) = (n — 1)!, 

T(x)^ V^expHeK-s, (38) 

dAbramowitz & Stegunll972h . which is valid for large x, x » 1 and 
gives sufficient accuracy (roughly 0.4% in the relevant Grange) for 
the purpose of this study. 




multipole order t 



Figure 4. Angular equilateral bispectra (y 3 ) (solid line), (ry 2 ) (dashed line) 
and (r 2 y) (dash-dotted line), as a function of multipole order t, for the red- 
shift disuibution of the DUNE galaxy sample. 
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Figure 5. Angular isosceles bispectra (y 3 ) (solid line), (ry 2 ) (dashed line) 
and (ry) (dash-dotted line), as a function of angle cost//, for multipole 
orders t = 10 2 (thick lines) and t = 10 3 (thin lines), for the redshift distri- 
bution of the DUNE galaxy sample. 



3.5 Angular bispectra 

This section gives an overview over the different configuration and 
scale dependences of the the angular bispectra (T q y 3 ~'>): The equi- 
lateral angular bispectrum B q {C, £, t) corresponding to the 3-point 
correlation functions {■f l y > ~ q ) are shown in Fig. [4] There is clearly 
a hierarchy in the bispectra, with <y 3 ) attaining the largest values, 
followed by (j 2 t) and (yr 2 ), separated by two orders of magnitude 
on large scales and up to seven orders of magnitude on small scales. 
The shape of the bispectra is determined by the mode-coupling 
functions M(k h kj) and the fc~ 2 -factors in perturbation theory. 

Apart from a clear scale dependence, the angular bispec- 
trum B q (i, cos iff) corresponding to the 3-point correlation func- 
tions (T q y 3 ~ q ) exhibit a different configuration dependence as well: 
Isosceles bispectra are shown in Fig. [5] as a function of the open- 
ing angle tf/ of the triangle. Typical variations of the bispectra with 



configuration on a fixed angular scale amount to two orders of mag- 
nitude. The configuration dependence itself changes with angular 
scale, as the bispectra plotted for £ = 100 and for £ = 1000 show 
different dependences on i/r. 

The configuration dependence of the bispectra (i J? y 3- ' 7 ) is 
shown for t 3 = 300 and £ 3 = 1000 i n Figs. [6] and [7] Specifi cally, I 
plot the configuration space variable dCoorav & S heth 2002) 



I 2 



B q {Ci,(2,£i) 



(39) 



which shows the dependence of the bispectrum B q on (fii, (%) with 
£1 fixed, in comparison to the equilateral configuration. The reason 
for this behaviour are the different power spectra of the density and 
potential fields, the latter having a much larger correlation length 
due to the AT 2 -factor. 



6 Bjorn Malte Schdfer 




multipole t\ 



Figure 6. Configuration dependence Rqfy(f\,£l) of the angular bispectra, 
(y 3 ) (bottom plane), (ry 2 ) (centre plane), (r 2 y) (top plane), as a function of 
t\ and (2 with (3 = 300, for the redshift distribution of the DUNE galaxy 
sample. The empty region to the left would contain invalid configurations 
violating the triangle inequality \£ 1 - t 2 \ < £3 < t\ + C 2 . 



£ = 10, 10 2 , 10 3 . The relative weakness of the non-Gaussianities 
in the iSW-effect is related to fact that r is dominated by contribu- 
tions from fluctuations on comparatively large spatial scales, where 
structure formation is well described by linear theory. 

In eqns. d40t and J41b . C TT (t) denotes the angular iSW power 
spectrum, C yy {t) the angular galaxy spectrum and C Ty {t) the cross- 
spectrum, respectively: 

J"XH 1 
<k — w 2 T (x)Q 2 M) P w (k = tlx), (42) 
X 

C^i) = (** d x 4 W y ( X )W T (x)Q y (a)Q T (a) P*(* = tlx), (43) 
Jo X 

Jr*XH 1 
dx - W 2 y (x)Q 2 y (a) P S s(k = £/ X ), (44) 
X 

with first-order time evolution Q y {a) and Q T {a) of the respec- 
tive fields (compare eqns. 1161 and|17|l and the weighting functions 
W T (x) and W y (x) (eqns. [I9]and|20j. P w (k) = P es (k)/(d„kf and 
P&p(k) = F 'ss(k) 1 '(d 2 H ) are the potential and the density-potential 
cross spectrum, where the Hubble distance d H = c/H makes the 
wave vector k dimensionless such that all three spectra have units 
of (Mpc/fc) 3 . 




multipole t\ 



multipole £ 2 



Figure 7. Configuration dependence Rg/ 3 (£i,£2) of the angular bispectra, 
(y 3 ) (bottom plane), (ry 2 ) (centre plane), (r 2 y) (top plane), as a function of 
l\ and I2 with £3 = 1000, for the redshift distribution of the DUNE galaxy 
sample. The empty region to the left would contain invalid configurations 
violating the triangle inequality \£ 1 - £2] ^ £3 < £\ + £l- 



3.6 Amplitude of the nonlinearities 

The overall amplitude of the nonlinearities probed by the iSW- 
effect can be quantified by the ratio a between the squared bispec- 
trum (e.g. for the equilateral configuration) and the corresponding 
power spectrum to the third power, 



a T {£) 



a y (£) = 



cut) 

B 2 ({, t,£)Y 
CU€) 



(40) 



(41) 



for which I obtain values of a T = 1.42x10^, 1.18xl0" 2 ,9.87xl0- 2 
and in comparison a y = 1.78 x 10 -1 , 1.98 x 10 -1 , 1.70 x 10" 1 for 



4 DETECTABILITY 
4.1 Covariances 

In contrast to power spectra, the observed bispectra B {a \, c) (l\,i2, 1-$) 
are unbiased estimates of the true bispectra B^ abc) {( l ,(2^i) m the 
case of Gaussian noise components, 



B(cibc)((l,i2,ti) — B( fl ) K ){tl,t2,(i). 



(45) 



Using a Gaussian approximation, the bispectrum covariance can be 
expanded into a product of three power spectra, and in the case of 
mixed bispectra, the covariance is a sum over all possible permu- 
tations of the observed cross-correlation C Ty (f ) and the two auto- 
correlations C TT (t) and Cyy(€). Specifically, the Gaussian approx- 
imation to the covariance CovfS^j^i,^, ^3), S(yt)(^,^,^)], 
which arises in the Wick-decomposition of the corresponding 6- 
point correlation function, assumes the shape (Hu 1999, 2002; 
iTakada & Jainll2003bl . l2004h 

Cov[S (fli)r) (A , t 2 , ti), B m {t\ ,f 2 J'i)] = 

C{ai)(£l)C(bj)({2)C{ck)(£l) 0(^5 1^6^ + 

C( a j)(ti)C(&£)(t2)C( CJ -)(t3) S[ l i' i 5e 2 i' } St 3 i' 2 + 

Qay)(A)C(io(^2)C((*)(^3) 5 1^6(^6^ + 

C(ai)(.£i)C(bi)(£2)C(cj)(h) 5t l e , 1 5e 1 t' l 8c 3 c[ + 

C(ak)(tl)C(bj)(l2)C(ck){tl) St^f Sl^St^ + 

C(ak){£i)C(bj)(£2)C(ci)(t?,) St^Si^Sf^ , (46) 

and scales like f^l for an observation which covers a fraction of 
/sky of the sky. C( fl ,)(Q corresponds to the observed CMB-spectrum 
C TT ({), if both indices are equal to one, the galaxy spectrum C yy (£) 
for both indices being equal to zero and the cross spectrum C Ty {€) 
for unequal indices: 



C C TT (£), a + i = 2, 
C (ai) (£) = \ C Ty (£), a + i=l, 
[ C yy (£), a + i = 0. 



(47) 
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N 


An 


/sky Z0 


b n 


3.0 x W 9 


2n 


0.5 0.64 


1.0 4.7 x10 s 



Table 1. Properties of the DUNE survey: total number N of objects, solid 
angle AO. covered (in radians), sky fraction /gw, redshift parameter zo, bias 
b and density per unit steradian n. 



I drop the negative sign of the cross spectrum C Ty (() because there 
will be always an even number of cross spectra in the terms of the 
expression for the covariance l |46t . 



4.2 Noise sources 

The observed spectra C TT (l) and C yy {i) differ from the theoretical 
expectations by the primary CMB fluctuations, a Gaussian instru- 
mental noise source <x; and a Gaussian beam B(C) in case of the 
CMB observation, and by a Poissonian noise term n in case of 
the galaxy survey, assuming that the noise sources are mutually un- 
corrected, 



C TT (0 = C TT (0 + C C MB(e) + a-;p-\£), 

C Ty (l) = C Ty {e\ 

Cyy(€) = Cyy({) + tl- [ . 



(48) 
(49) 
(50) 



More specifically, the observational noise consists of these contri- 
butions: 

(i) The Fourier-transform of a Gaussian beam is given by 
p- 2 (t) = exp(A6 2 £(£ + 1)). For the beam width I use the value 
A9 = T,\ corresponding to the v = 143 GHz channels closest to 
the CMB emission maximum. The value Tcmb - 2.725 K for the 



CMB temperature is used when converting w 



3 cri to the 



noise amplitude in the dimensionless temperature perturbation r, 
with w- 1 = (0.01 |i.K) 2 fcaldarriagaetal.il 19971) . 

(ii) In addition, I generated a CMB temperature power spec- 
trum CcmbM, equally scaled with the CMB temperature 7cmb = 
2.725 K, with the Code for Anisotro pies in the Microwave Back- 
ground (CAMB, iLewis et alj[2OO0t) for the fiducial ACDM cos- 
mology. The noise contribution from the CMB-spectrum CcmbOT 
is the main difficulty in observing the iSW-bispectra, because it pro- 
vides high values for the covariance at low multipoles £, and it by 
far dominates C TT (f), Ccmb(^) S> C tt (£) on the angular scales con- 
sidered. 

(iii) The Poissonian noise term in the galaxy counts is the in- 
verse of the number density n of objects per unit steradian. TableQ] 
summarises properties of the DUNE main galaxy sample. The main 
advantage of DUNE is the large sky coverage and the high num- 
ber of objects. For simplicity, I have assumed a constant (i.e. non- 
evolving) unit bias for the DUNE galaxy sample. 

An overview over the observed spectra is given by Fig. |H in com- 
parison to the noiseless spectra: For the angular scales considered 
here, the Poisson noise term in the galaxy number counts is not yet 
an issue, the cross-spectrum and the intrinsic CMB spectrum are of 
similar magnitude, and at £ = 10 3 about a factor of 10 6 larger than 
the noiseless iSW-spectrum. 



|http : //c amb . info/ 
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Figure 8. The observed spectra C y y(t) (thick solid line), C Ty (l) (thick 
dashed line) and C TT (C) (thick dash-dotted line), for PLANCK in conjunc- 
tion with DUNE. In comparison, the corresponding noiseless spectra C yy (€) 
(thin solid line) and C TT (f) (thin dash-dotted line) are also plotted. 



4.3 Signal to noise ratios 

The signal to noise ratio S ? (assuming a Gaussian likel ihood) for 
a measurement of the bispectrum (r' 7 y 3 ~ ? ) is g iven by jHull200rA 
iTakada & Jainll2003bl . |2004 iDolnev et alj|200fj) : 

trm bm 

^= Z Z Z Z w^)(&p.4)%^p- ( 51 ) 

lp=tmn t'=e min (*) SM 
1 a+b+c-q ('+ j-^k-q 

with ( p e {(1,(21^3} an d t'p 6 {^p^2'^}' F° r conven i ence > I abbre- 
viated the inverse covariance of the bispectrum B( a &c)(^i>^2,^3), 

G2#p. 9 = (Cov[B (flfc) (A , £ 2 , £ 3 ), B m (€\ , £' 2 , . (52) 

Like lTakada & Jainl d2003bll2004l) . I use binned summations in the 
multipoles t\ and l 2 , but carry out an unbinned summation in £ 3 
in order to account for the vanishing Wigner-3 j symbol if Yi P ip 
is an odd number and for the sign change of the Wigner-3 j symbol 
depending on whether t p modulo 4 vanishes or not. Fig.|9]shows 
the cumulative signal to noise ratio of the measurement of B q as a 
function of maximum multipole order f max . In the following I will 
use kt\ = Mi = 30 for the largest £ max -values considered, and the 
summation is carried out starting from f m j„ = 10. When computing 
the covariance of the bispectrum (T q y 3 ^' 1 ), one needs to make sure 
that each term in the summation contains exactly 2q iSW-fields r, 
which is achieved by restricting the summation over all possible 
permutations to those terms which satisfy the conditions a + b + c = 
q and ( + j + k = q. 

As shown by Fig. [9] the cumulative signal to noise ratio ~L q 
increases slowly with increasing multipole order £. The galaxy bis- 
pectrum (y 3 } should be measurable with 3cr significance on de- 
gree angular scales. The bispectrum (ry 2 ), sadly, only reaches a 
significance level of 0.4<r up to the scale £ = 10 3 , and one 
would need to carry on the measurement to angular scales of =s 10 4 
for a detection significance of 3cr. These angular scales are be- 
yond PLANCK's resolution limit and other secondary anisotropies 
such as the Sunyaev-Zel'dovich effect or the contribution to point 
sources becomes important. Furthermore, the second order pertur- 
bation used in this work would be no longer applicable. At £ = 10 3 , 
the So is roughly two orders of magnitude larger than Si, which in 
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10' 

multipole € 

Figure 9. Signal to noise ratio £„ of the bispectrum (r , y 3_<? ) as a func- 
tion of maximum multipole moment f m - dx : Xo (solid line), Ej (dashed line) 
E2 (dash-dotted line), for cross-correlating the DUNE galaxy sample with 
PLANCK CMB data. The horizontal lines mark the confidence levels lcr 
(thick dashed line) and 10cr (thick solid line). 



<r 3 > 



(ry 2 ) 



(r 2 y) 



So = 56.29 



: 0.593 £ 2 = 1.667xl0- 



Table 2. Cumulative signal to noise ratios Z 9 for measurements of the three- 
point correlation functions (^y 3 ' 9 ), q = 0, 1, 2, for PLANCK CMB data in 
cross-correlation with the DUNE survey, up to PLANCK's resolution limit 
&™ = 2x 10 3 . 



turn is slightly more than three orders of magnitude larger than Z 2 . 
S 2 is separated from £3 by four orders of magnitude. This behaviour 
reflects the increasingly higher powers of the primary CMB power 
spectrum Ccmb(^) which enters the covariance as a noise source. 

Signal to noise ratios ~L q for measurements of all three-point 
functions by correlation with the DUNE main galaxy sample are 
compiled in Table [2] where the summation was extended to the 
maximum multipole moment i = 2x 10 3 . The value of the signal to 
noise ratio Si for a measurement of the bispectrum (ry 2 ) amounts 
to 0.6<x considering triangle configurations up to PLANCK's reso- 
lution limit, which is not large enough to be realistically detected, 
because the confidence level would only be erf (0.6/ V2) - 0.45. As 
a consistency check, I also derived the signal to noise ratio £3 for 
measuring the bispectrum (r 3 ), and obtain a value of £3 - 10~ 8 cr, 
which is quite comparable to the estimate of lCooravl (|2002), keep- 
ing in mind that we use different cosmological models and a differ- 
ent perturbative approach. 



4.4 Degeneracies 

The Fisher matrix jTegmark etaill 19971) . which describes the log- 
arithmic decrease in likelihood if the cosmological parameters 
move away from their fiducial values, is defined in analogy to 
eqn. {5J}: 



{abc) {ijk) 
a+b+c-q i+j+k-q 



dB (abc) (£ p ) ,. k 



dB Wk) ({' p ) 



(53) 



but the derivation of parameter degeneracies would not yield com- 
petitive parameter bounds, given the low values for the signal to 
noise ratio ~L q . Quite generally, all bispectra are proportional to cr^ 
because of the perturbative evaluation to second order, and oc Q? m 
because the matter density plays the role of a coupling strength in 
the line of sight expression for the iSW perturbation r. The bis- 
pectra with small q are more strongly influenced by the Hubble 
parameter h via the shape parameter T = Cl„,h because the strong 
weighting oc 1/fc 4 of the potential power spectrum dampens the 
sensitivity. Concerning the line of sight integrated quantities the 
scale factor a at which equality Q.„,(a) = ClA(a) is reached, plays 
an important role, because at the corresponding redshift most of 
the iSW-signal is created. Furthermore, models with small average 
dark energy equation of state parameter w give rise to an iSW-effect 
at higher redshifts compared to ACDM. 



5 SUMMARY 

The aim of this study is an investigation of the detectability of the 
nonlinear iSW-effect using mixed bispectra of the form <T ? y 3_? ), 
q = 0, 1, 2, between the galaxy density 7 and the iSW-temperature 
perturbation r. The bispectra were consistently derived in second 
order perturbation theory. I investigated the time evolution of the 
souce terms, the configuration dependence of the bispectra, the 
achievable signal to noise ratio in a measurement cross-correlating 
PLANCK and DUNE data, and their parameter sensitivity. 

(i) The nonlinear iSW-effect has a different parameter depen- 
dence compared to the linear iSW-effect, as it is sensitive to the 
derivative d(.D 2 /a)/da instead of d(D + /a)/da. Particularly, it does 
not vanish in SCDM-models, where Q.„, = 1 and consequently 
D + (a) = a. 

(ii) I employed second order hyper-extended perturbation the- 
ory for deriving the the bispectra, and used a Gaussian approxima- 
tion for describing the covariance of the measurements. As noise 
sources, I considered the intrinsic CMB fluctuations, pixel noise 
and a Gaussian beam for the CMB observation, and a Poissonian 
noise term in the galaxy density. To this point, I worked with a 
constant linear biasing model for relating the fluctuations in galaxy 
number density to those of the dark matter density. 

(iii) The configuration and scale dependence of the mixed bis- 
pectra reflects the interplay between the correlation length of the 
density field and the much larger correlation length of the gravita- 
tional potential. 

(iv) I computed the signal to noise ratio for the measurements 
of bispectra (r^y 3-9 ), q = 0, 1, 2, with a Gaussian approximation to 
the covariance, in which the intrinsic CMB fluctuatios are the most 
important noise source which make the bispectra difficult to ob- 
serve. Values for the cumulative the signal to noise ratio = 0.6cr 
for the measurement of the bispectrum (ry 2 ) are obtained for cross 
correlating PLANCK data with the DUNE main galaxy sample, up 
to PLANCK's resolution limit at I = 2x 10 3 , where the dominating 
noise source are the primordial CMB fluctuations. 

(v) An algorithm for evaluating 3-point correlation function 
would need to be implemented for PLANCK data processing. This 
is most likely an algorithm operating in harmonic space instead of 
real space, because the evaluation of the 3-point correlation func- 
tions scales as A^ 3 ix where N vix is the number of pixels, which be- 
comes prohibitive at N pix ^ 10 7 , corresponding to multipole orders 
of I = 1000. From the observational point of view, the influence 
of unresolved microwave point sources or the kinetic Sunyaev- 
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Zel'dovich effect, which are equally associated with overdensities 
in the LSS, on a measurement of the iSW- temperature fluctuation 
is yet unquantified. 

Future studies will treat three-point correlation functions of 
the type {f 1 ^' 1 ) between the iSW temperature perturbation r and 
weak gravitational lensing convergence k, which would be attrac- 
tive because lensing measures directly the fluctuations in the dark 
matter density, without uncertainties related to bias and bias evolu- 
tion. In comparison to (T c 'y i ^ q ), the bispectrum (i^k 3 ^) would mea- 
sure the dark energy properties at lower redshifts. For the DUNE 
sample, the constraints would come from redshifts of z — 0.4, com- 
pared to z - 0.9, with more strongly evolved nonlinear structures. 
A third point in favour of gravitational lensing is the smaller sam- 
pling noise in the galaxy ellipticity cr\ln - 0.1 jn compared to 1/n 
in the galaxy density. 
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